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£^ | Abstract. We have evaluated forward dispersion relations for scattering amplitudes that follow 

| • from direct fits to several sets of 7tn scattering experiments, together with the precise K decay 

' results, and high to energy data. We find that some of the most commonly used experimental sets, as 

well as some recent theoretical analyses based on Roy equations, do not satisfy these constraints by 
— ■ several standard deviations. Finally, we provide a consistent n% amplitude by improving a global fit 

' to data with these dispersion relations. 
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CS ; INTRODUCTION 

A precise knowledge of the %% scattering amplitude has become increasingly important 
since it provides crucial tests for one and two loop Chiral Perturbation Theory (ChPT), as 
well as crucial information on three topics under intensive experimental and theoretical 
investigation: light meson spectroscopy, pionic atom decays and CP violation in kaons. 
Unfortunately, these precision studies are very cumbersome due to the poor quality of the 
data which is affected by large systematic errors. Here we review our recent work where 
we checked the fulfillment of dispersion relations by different sets of data commonly 
used in the literature, and provided parametrizations consistent with such requirements. 

Recently, Ananthanarayan, Colangelo, Gasser and Leutwyler (ACGL)ll|] and Colan- 
gelo, Gasser and Leutwyler (CGL)[2] have used data, analyticity and unitarity through 
Roy equations, and ChPT, to build a %% amplitude. They provide phase shifts up to 
0.8 GeV, scattering lengths and effective ranges claiming an outstanding precision. 
While the methods of CGL constitute a substantial improvement over previous ones, 
their analysis has to rely on some input, part of which we have recently questioned. First 
of all, their Regge high energy representation does not describe the high energy data, 
and does not satisfy well certain sum rules. Second, their D2 wave is incompatible with 
a number of requirements. Finally some of their input has remarkably small errors and 
relies precisely on some data sets that do not satisfy well forward dispersion relations. 
All this is discussed in the present note, which is based our recent works J3l0,0]. 
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FIGURE 1. a) P wave phase shifts from TtTt scattering sua compared with the prediction with the 
parameters (2.1) (solid line below 1 GeV). Note that this is not a fit to these data, but is obtained from 
the pion form factor |6|. The error here is like the thickness of the line. Above 1 GeV, the dotted line and 
error (PY) are as follows from the fit in 131. b) SO phase shifts and error band as given by Eqs.( I5I6> below 
1 GeV, and from [5] above. The K/4 and K% % decay data are not shown, (see our |5] for details). 

LOW ENERGY PARTIAL WAVES FROM FITS TO DATA 

The SO, S2 and P partial waves at low energy, s 1 / 2 ^ 1 GeV 

We first consider wave-by- wave fits to data for the SO, S2, P waves, as in o], which 
improve our "tentative solution" in ||3|]. To fit the phase shifts, 8 (s), we parametrize 
cot 8 (s) taking into account its analytic properties, as well as its zeros (associated with 
resonances) and poles (when the phase shift crosses tin, n =integer). 

For the P wave, up to ~ 1 GeV we use the results from a fit to the pion form factor as 
given in The comparison with %% scattering data can be seen in Fig.l. We take so as 
the point at which inelasticity begins to be nonegligible, and we write 

cot 8 1 (s) = ^ (Mj - s) ji?o + B ! ^ ~ ^Z! | ; s l J 2 = 1.05 GeV. (1) 

B = 1.069 ±0.011, Bi =0.13 ±0.05, M p = 773.6±0.9, 

ai = (37.6 ±1.1) x lO^M^ 3 , b\ = (4.73 ±0.26) x lO^M^ 5 ; (2) 

For the S2 wave we fit data where two like charge pions are produced: [7] although 
these pions are not all on their mass shell, at least there is no problem of interference 
among various isospin states. At low energies, we fix the Adler zero at zi = M n and fit 
only the low energy data, s 1 / 2 < 1.0 GeV; later on we allow 12 to vary. We have 

cot^'M = * U +Bl f 7 f^ZZ } , «,'/* = 1.05 OeV, (3) 
B = -80.4 ±2.8, Bi = -73.6 ±12.6; 



a { Q 2) = (-0.052 ±0.012)M~ 1 ; bf 1 = (-0.085 ±0.01 1)M^ 3 . (4) 

The SO wave experimental situation is somewhat confusing, and we consider two 
methods of data selection. In both, we fit K14 and K — > 2k decay data, in which pions 
are on the mass shell. In the first method, called global fit, we include some points at 
0.81 GeV < s 1 / 2 < 0.97 GeV, where the various experiments agree within ^ 1.5 a. 
Care is exercised to compose errors realistically, see details in [5], Subsect 2.2.2. In this 
case we fix the Adler zero at zo = M K and find 

•(0)/^_s 1/2 M l M l- s f D , D \/so--s 



B = 21.04, =6.62, M a = 782±24MeV; 8^°\m K ) = 41.0° ±2.1° ; 
aj) 0) = (0.230 ±0m0)M K \ b { ° ] = (0.268 ± 0.01 1)M K 3 ; 

this fit (shown in Fig l.b as PY) is valid for s 1 / 2 < 0.95 GeV. The Bj errors are strongly 
correlated; uncorrelated errors are obtained if using the parameters x, y with 

B =y-x; fli = 6.62-2.59*; y = 21.04±0.70, x = 0±2.6. (6) 

The other method is to fit only K14 and K — ► 2n data, or to add to this, individually, 
data from the various experimental analyses. The results can be found in Table 1. 



The SO, S2 and P partial waves at 1 GeV <l s 1 / 2 <^ 1.42 GeV 

The D and F data are scanty, and have large errors. To stabilize the fits we impose 
the values of the scattering lengths that follow from the Froissart-Gribov representation. 
This is not circular reasoning since their Froissart-Gribov representation depends mostly 
on the SO, S2 and P waves, and very little on the DO, D2, F waves themselves. We do 
not discuss here the DO and F waves (see [5]) as they do not present special features. 

For D2 we only expect important inelasticity when the nit — > pp channel opens up, 

so that so = 1.45 2 GeV 2 ~ 4M 2 ,. A pole term is necessary here, since we expect 8^ to 

1— I (2) 

change sign near threshold: the data [7] give negative and small values for 8 2 above 



some 500 MeV, while, from the Froissart-Gribov representation, it is known[ll] that 
the scattering length must be positive. Indeed we include in the fit the value a\ = 
(2.72 ±0.36) x 10 _4 M^ 5 . In addition, the clear inflection seen in data around 1 GeV 
asks for a third order conformal expansion. So we write 

cot8 2 (s) = w {B Q+ BMs)+B 2 w(s) ) - - - , w(s) = ^ + 

And we find B = (2.4 ± 0.3) x 10 3 , B x = (7.8 ± 0.8) x 10 3 , B 2 = (23.7 ± 3.8) x 10 3 , 
A = 196 ± 20 MeV. The fit, which may be found in Fig 2, returns reasonable numbers 

for the scattering length and for the effective range parameter, b 2 ■ 

4 2) = (2.5±0.9) x 10" 4 M^ 5 ; bf ] = (-2.7 ±0.8) x 10- 4 M,r 7 . (7) 
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FIGURE 2. Continuous line: The I ~2, D-wave phase shift, obtained by only fitting the experimental 
data. Broken line: with the parameters improved using dispersion relations. Dotted line: the fit, valid 
between i 1 / 2 =0.625 GeV and 1.375 GeV, of Martin, Morgan and Shaw which ACGL and CGL, however, 
use from threshold to s l l 2 = 2 GeV. The experimental points are from |7]. 



THE HIGH ENERGY (S 1 / 2 > 1.42 GEV) INPUT 

In order to test dispersion relations we also need the imaginary part of the scattering 
amplitude at s 1 / 2 > 1.42 GeV, that we take from a Regge fit to data [4] (and the slightly 
improved rho residue of [J]). We note that, in the early 1970s, when nn phase shifts 
were poorly known and, above all, when it was still not clear that the standard Regge 
picture is a QCD feature, Regge factorization was questioned ifn using crossing sum 
rules and then-existing low energy phase shift data. This was adopted later by ACGL and 
CGL, assuming a too large rho residue and a Pomeron a third of what factorization and 
the experimental data on the total ktz cross section implies, as well as unconventional 
slopes. Unfortunately this has been also used in subsequent Roy equation analyses. As 
discussed in 0,0], however, standard Regge factorization describes experiment ifLlLfLlll 
and is perfectly consistent with crossing sum rules if assumed to hold above 1.42 GeV. 

In Fig 3 we show our Regge description of the imaginary parts of nn scattering 
amplitudes USUI] together with the dataJH], compared with that used by ACGLfl], 
CGL|H above 2 GeV. 

Between 1.42 GeV < s 1 / 2 < 2 GeV, these authors use the scattering amplitude re- 
constructed from one Cern-Munich phase shift analysis, and, in particular for SO, the 
re-elaboration of Au, Morgan and Pennington II 1(11 . Unfortunately, in this region the 
inelasticity is large and the Cern-Munich experiments, which only measure the differ- 
ential cross section for %% — > %% are insufficient to reconstruct without ambiguity the 
full imaginary part. In addition, in ^ Ifl] we showed that Cern-Munich phases fail to 
pass a number of consistency tests. This is also seen clearly in Fig 3, where we plot the 
total cross section for that follows Hyams et al.- llOll - which is incompatible with 

other experimental data as well as with Regge factorization. 
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FIGURE 3. The %% cross sections. Experimental points from 1 13]. The stars at 1.38 and 1.42 GeV (PY) 
are from the phase shift analysis of experimental data given in |5]. Continuous lines, from 1.42 GeV (PY): 
Regge formula, with parameters as in jj] (the three lines per fit cover the error in the theoretical values 
of the Regge residues). Dashed lines, above 2 GeV: the cross sections following from ACGL;| 1 ] the gray 
band covers their error band. Below 2 GeV, the dotted line corresponds to the n + n~ cross section from 
the Cern-Munich analysis; cf. Fig 7 in the paper of Hyams et al. [ 10] 



In a recent paper Caprini, Colangelo, Gasser and Leutwyler. lfl7ll to be denoted by 
CCGL, review our work in Jj] and conclude that, still, they consider the CGL solution 
consistent. They also raised the contention that our Reggeistics could not be correct 
because it violates certain sum rules. In view of Fig. 3 this contention is meaningless 
since the PY cross sections are perfectly compatible with high energy (s 1 / 2 > 1.42 GeV) 
experimental data, while the ACGL ones are not. In O, |j, |5|] we also checked that our 
representation satisfies two crossing sum rules. 

Concerning D2, ACGL and CGL borrow an old fit in the book of Martin, Morgan and 
Shaw. lfl5ll where only intermediate energy data were fitted. 

S 2 (2) (s) = -0.003(5/4M 2 ) (1 -AMl/s) 5 '\ (8) 

which fails at threshold (it gives a negative scattering length) and does not fit well data 
below 1.42 GeV, as shown in Fig 2. Above 1 GeV, this D2 phase grows quadratically 
with the energy, while Regge theory predicts all phases to go to a multiple of n. In 
particular D2 should go to zero; see Appendix C of yD for details. It is true that this D2 
wave is small but, given the accuracy claimed by CGL, it is certainly not negligible. 



CHECKING FORWARD DISPERSION RELATIONS 



In the present Section we study how well the previous amplitudes obtained from fits to 
different sets of data satisfy forward dispersion relations. We consider three indepen- 
dent scattering amplitudes in ^-symmetric or antisymmetric combinations, that form a 
complete set: — > — > and the t channel isospin one amplitude, 

I t = I. The reason is that the two first depend only on two isospin states, and have posi- 
tivity properties: their imaginary parts are sums of positive terms, thus reducing the final 
uncertainties. Hence, for we have 

ReFoo(s) -F Q0 (4M n ) —P.P. J^ds s , {s , _ s)(s , _ 4M 2 )(s , + s _ 4M 2y P) 

In particular, for s = 2M%, which will be important for the Adler zeros, we have 

F 00 (4Ml) = F 00 (2Ml) +D 00 , A)0 = ^ f 2 ds (10) 

% JamI s{s-2M£){s-4M£) 

For the n°7t + channel, which does not depend on SO: 

RcF 0+ (s)-F 0+ (4M n ) P.P.y^d, sW _ AM 2 ){sl + s _ AM 2y 

At the point s = 2M\, this becomes 

Finally, for isospin unit exchange, which does not require subtractions, 

{ ' J k hul (s'-s)(s' + s-4MlY {U) 

at threshold this is known as the Olsson sum rule. 

Depending on the method we use to fit the SO wave we find the results in Table 1, 
where, we have separated on top those fits to data with a total x 2 /d-o.f. < 6 for the 
7r°7r° and I t = 1 dispersion relations up to 0.925 GeV, a fairly reasonable j} jd.o.f. 
since these fits were obtained independently of the dispersive approach. 

However, in Table 1 we also list the very frequently used t and s-channel solutions 
of Estabrooks and Martin [10], those of Protopopescu et al. [9], from Table VI, VIII and 
table XII, as well as the solution A of Grayer et al. [ 10]. Their I t = 1 plus dispersion 
relation total % 2 /d.o.f. is surprisingly poor: 11.3, 10.1, 7, 6, 7.5, 9.9, respectively. 
Therefore, any result that relies heavily on these sets should be taken very cautiously. 



TABLE 1. PY: our global fit, Eas. J5l6> . We do not give its Bo an d B\ uncertainties as they are 
strongly correlated, see Eq.© for the uncorrelated ones. Grayer B, C, E: different solutions in Grayer 
et al.CB- Kaminski: 0U- InH we have also studied fits to the data in Tables VI, XII and VIII in 01, 
to Solution A in[ 10], as well as fits to the theoretical outcome in Estabrooks and Martin.] 10]. They all 
give a total %' d.o.f. > 6 
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IMPROVED FITS USING DISPERSION RELATIONS 

We now improve the previous low energy fits parameters by fitting also the dispersion 
relations up to 0.925 GeV, thus obtaining parametrizations more compatible with ana- 
lyticity and s — u crossing. This is an alternative method to Roy equations; it is better in 
that we do not need as input the scattering amplitude for \t \ up to 30M 2 ., where the Regge 
fits existing in the literature disagree strongly (see [5], Appendix B) and also in that we 
can test all energieso], whereas Roy equations are valid for s 1 / 2 < y/60M K ~ 1.1 GeV 
(and only applied up to 0.8 GeV). Starting from Eqs. (l5l6l) . we find, in M n units, 

SO; s 1/2 < 2m K 



S2; s l/1 < 1.0 



B = 17.4±0.5; B x =4.3 ±1.4; 
M(j = 790 ±21 MeV; zo = 195 MeV [Fixed]; 
4 0) =0.230 ±0.015; = 0.312±0.014. 
B = -80.8 ±1.7; B x = -77 ±5; z 2 = 147 MeV [Fixed]; 
4 2) = -0.0480 ±0.0046; b { 2) = -0.090 ±0.006. 



S2; 1.0<s 1/2 < 1.42: B = -125 ±6; B 1 = -119±14; e = 0.17±0.12. 

P; s l/2 < 1.05: B = 1.064±0.11; B x = 0.170 ±0.040; M p = 773.6 ±0.9 MeV; 
a\ = (38.7 ±1.0) x 10~ 3 ; b x = (4.55 ±0.21) x 10~ 3 . 
DO; s l/2 < 1.42: B = 23.5 ±0.7; B x = 24.8 ± 1.0; £ = 0.262 ±0.030; 

4 0) = (18.4 ±3.0) x 10~ 4 ; bf ] = (-8.6 ±3.4) x 10~ 4 . 
D2; s l/2 < 1.42: B = (2.9 ±0.2) x 10 3 ; B x = (7.3 ±0.8) x 10 3 ; 

B 2 = (25.4±3.6) x 10 3 ; A = 212 ± 19; 
4 2) = (2.4 ±0.7) x 10~ 4 ; b {2) = (-2.5 ±0.6) x 10~ 4 . 
F; s l/2 < 1.42 : B = (1.09± 0.03) x 10 5 ; B x = (1.41 ±0.04) x 10 5 ; 

a 3 = (7.0±0.8) x 10" 5 . (13) 

In Fig. 4 we show the improved curves for SO and S2, and that of D2 in Fig. 2. 




FIGURE 4. a) The improved SO phase shift (PY improved, Eq|l3}, the global fit (PY from data, he SO 
Eqs. J5I6» . and the improved solutions "K decay only" and "Grayer C" of Table 2 (almost on top of PY 
improved). The solution CGL||2|] (dashed line) is also shown, b) S2 improved Phase shift (PY improved, 
Eq. global fit (PY from data, Eq. @); the solution CGL yj (thin continuous line) and the improved 
parametrization with K decays and So. B of Grayer et al. fToll . 

Concerning the improved fits to individual sets of data, we get somewhat different 
results for SO, listed in Table 2. In Table 2 we also show the % 2 /d.o.f. of each forward 
dispersion relation and the standard deviations for the sum rule in EqUQl which are 
more than four for K decay plus the Grayer B or E or Kaminski improved solutions. 
Concerning the other waves, no matter what set of parameters from data fits we start 
from, we end up with very similar values to those given in Eq{T21 This can be checked 
in Fig.4.b, where we show the improved "K decay + Grayer Sol. B" S2 wave. Even 
though it is the one for which we obtained the most different central values for the SO 
wave compared with those given in Eq{T3j it falls perfectly within the uncertainty of our 
improved solution. 



TABLE 2. Improved fits. Names are as in Table 1 . Although errors are given for the Adler zero, we fix 
it when evaluating other errors, to break the otherwise very large correlations 
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DISPERSION RELATIONS AND THE CGL SOLUTION 

We have also checked the fulfillment of forward dispersion relations for the CGL 
solution for the SO, S2 and P waves al low energy. This is depicted in Fig 5, where 
we show, both for CGL and our improved fit, Eq.dTjT). the mismatch between the real 
part and the dispersive evaluations, that is to say, the differences A,-, 



Ai=ReF^(j,0) -P.P. / ds'— — — v ' (14) 

v J % J ml (s'-s){s'+s-AMiy 



Aoo = ReFoo(5)-F o(4M^) (15) 

_s(s-4M^ [°° {2s'-AMl)lmF m {s') 

K ' " J AMI -s)(s> -4Ml)(s> + s -4MD' 

A 0+ = RcF 0+ (s)-F 0+ (4Ml) (16) 

4s-*Mp r» {2s'-AMl)lmF 0+ (s') 

K ' 'J AMI s'(s'-s)(s'-4Ml)( S ' + s-4Ml)' 

These quantities would vanish, A,; = 0, if the dispersion relations were exactly satisfied. 

We include in the comparison of Fig 5 the uncertainties; in the case of CGL, 
these errors are as follow from the parametrizations given by these authors in |2y, for 
s 1 / 2 0.8 GeV. At higher energies they are taken from data via our parametrizations. 
It can be clearly seen that the CGL parametrizations do not satisfy these forward disper- 
sion relations by several standard deviations. 

One might wonder why in Table 2, the SO wave improved "K decay+solution B" 
yields % 2 /d.o.f. of order one for the forward dispersion relations, being so similar to 
what CGL get for that wave. The reason is that, as we show in Fig.4.b. the improved 
solution B, requires an S2 wave that is even farther from the CGL S2 wave than our 
global improved solution. 
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FIGURE 5. Dispersion relations for the 7in amplitudes of |2] (CGL) and for our improved global fit 
(PY improved, Eq^] We plot the differences A,-, Eqs.^J between real parts calculated directly from the 
parametrizations, or from the dispersive formulas. Consistency within one sigma occurs within the shaded 
bands. The progressive deterioration of the CGL results as the energy increases is apparent here. 



LOW ENERGY PARAMETERS IN THE LITERATURE 

We here present, in Table 3, the low energy parameters obtained from Roy equations by 
CGL, by Descotes et al Jl5ll . that we denote by DFGS, and by Kamihski, Lesniak and 
Loiseau lll5ll . denoted by KLL. This is compared with what we find fitting experimental 
data, improved with dispersion relations (see [5] for details). 

The mismatches between many of the parameters of CGL and PY are apparent here; 
not surprisingly, they affect mostly parameters sensitive to high energy: 

• The S wave parameters are compatible, mainly due to the large PY uncertainties. 

• The mismatch between CGL and PY for and is roughly 2.5c. In ref.j3|] 
we pointed out that the ACGL and CGL results did not satisfy the Froissart Gribov 
sum rules. This happens to more than four standard deviations for the difference 
between the CGL calculation using Wanders sum rules minus the Froissart-Gribov 
representation. This larger mismatch, as pointed out in ifTvh - does not involve the S 
and P waves, and is due to the Regge and L > 2 wave input, evidently very different 

from the beginning for CGL and PY, but that certainly affects the values of ai an d 

• The b\ calculation differs by more than 4 standard deviations for PY and CGL. 



TABLE 3. Units of M % . The numbers in the CGL column are as given by CGL in Table 2 and 
elsewhere in their text. In PY, the values for the D, F waves parameters are from the Froissart- 
Gribov representation. The rest are from the fits, improved with dispersion relations, except for a\ 
and b\ that have been taken as in |5[. 
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THE ACGL, CGL PHASE AT S 1 / 2 = 0.8 GEV 

In the ACGL, CGL analyses, by the input phases for the SO, S2 and P waves at the point, 
s 1 / 2 = 0.8 GeV, where they match the solutions to the Roy equations to the experimental 
amplitude. Indeed it is dominant for their Olsson sum rule calculation, which involves 
the I t = 1 channel. 

The quantity 5^ 0) ((0.8 GeV) 2 ) is in fact given in Eq. (7.3) of ACGL as 

5,S 0) ((0.8 GeV) 2 ) =82.3 ±3.4°. (17) 

whose error may be contrasted with the estimates of |@], which vary, for the data above 

0.8 GeV, between 6° and 18°, or with the 5^ 0) ((0.8 GeV) 2 ) values we obtain from fits 
to different sets of data in Table 1, or the improved fits in Table 2. Small errors could 
be expected from theoretical analysis including many data but the above small error was 
used as an input. 

The reason to consider such an small error is that ACGL consider the difference 
8\ — 8q ^ at 0.8 GeV, in the hope that some of the uncertainties will cancel. Then 
they interpolate and then average the points from a choice of three analysis of the 
CERN/Munich experimentlfioll 

r 23.4 ±4.0° [Hyamsetal.] 
Si ((0.8 GeV) 2 ) - 5 (0) ((0.8 GeV) 2 ) = < 24.8 ± 3.8° [Estabrooks and Martin, s-channel] 

I 30.3 ± 3.4° [Estabrooks and Martin, r-channel]. 

and set 

5 (0) ((0.8 GeV) 2 ) - Si ((0.8 GeV) 2 ) = 26.6 ±2.8°. 



However, this error does not include systematics. All numbers here stem from the same 
experiment, and differ only on the method of analysis. Their spread is an indication 
of the systematic uncertainties, roughly an additional ±4° . In addition, the Hyams et 
al. value above is only one of five solutions in Grayer et al. and considering also 
data of Protopopescu et al.,[9] the systematic error would increase to 10° . Remarkably, 
Estabrooks and Martin themselves, point out in their section 4 (first paragraph) that 
different D wave input "lead to systematic changes in 8® of the order of 10°". 
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